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We present our recent results for the B-parameters, which parameterize the AB=2 transition amplitudes. 
Calculations are made in quenched QCD at /9=5.7, 5.9, and 6.1, using NRQCD for heavy quark and the 
0(a)-improved action for light quark. The operators are perturbatively renormalized including corrections of 
0{a s /amo). We examine scaling behavior of the _B-parameters in detail, and discuss the systematic uncertainties 
using scatter of results with different analysis procedures adopted. As a result, we find BB d (rrib) = 0.84(2) (8), 
B B JB Bd = 1.017(10)(±o) and Bs 3 (m b ) = 0.87(1) (9) in the quenched approximation. 



1. Introduction 

The precise determination of AB=2 transi- 
tion amplitudes would allow us to explore new 
physics beyond the Standard Model. A number 
of works to compute the amplitudes have been 
already done and are still in progress using the 
lattice QCD ||. The AB=2 transition ampli- 
tudes are conventionally expressed in terms of 
the corresponding B-parameters, which param- 
eterize the deviation from the vacuum satura- 
tion approximation of the amplitudes. Here we 
present our recent results on the two phenomeno- 
logically important B-parameters, Bs q and Bs q , 
which are defined in the AB=2 transitions with 
O l = huPLq h^Lq and O s = bP L q bP L q, re- 
spectively, where q=s or d and Pl = 1 — 75. 

One of the features in this work is the use of lat- 
tice NRQCD H for heavy quark, which enables 
us to obtain a better control of a systematic error 
than other approaches since it does not need any 
extrapolation in heavy quark mass. The price 
one has to pay, however, is that in using effec- 
tive theory one can not take the continuum limit 



and so that the power-law divergence of the form 
a™/(amb) 1 and the discretization error of o 2 Aq CD 
appear simultaneously. This means that we have 
to choose lattice spacing such that both sources 
of systematic errors are under control. To see this 
we perform quenched QCD simulations at [3=5.7, 
5.9, and 6.1 and investigate stability of the results 
against the lattice spacing. 

The simulation parameters are summarized in 
Table |. We use the standard plaquette action for 
gauge field and the 0(a)-improved action for light 
quark. In order to see the chiral behavior, we take 
four different masses for light quarks in the range 
[m s /2, m s ). We also take five different masses for 
heavy quarks in the range [2m&/3, 4m;,] to inves- 
tigate the convergence of 1/wq expansion. The 
lattice spacing for each lattice, shown in Tabic [l], 
is determined by the string tension. 

2. Method 

The operators Ox in the continuum are written 
in terms of those on the lattice Oj^ as 
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6.1 


5.9 


5.7 


size 


24 a x 64 


16 3 x 48 


12 3 x 32 


# of conf. 


518 


419 


420 


l/a c [GeV] 


2.29 


1.64 


1.08 



Table 1 

Lattice parameters in this work. 

The one-loop matching coefficients zx.y (X=L, 
S) as well as the precise definitions of operators 
and other details can be found in Ref . || , where 
the perturbative calculations are performed in- 
cluding 0(a s / amo) corrections. 

The b quark field defining each operator such 
as Ol and Os is related to the non-relativistic 
field h—(Q x T through the inverse of Foldy- 
Wouthuysen-Tani transformation as 
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It should be noted that because of a lack of 
the matching coefficients at the higher order 
of 0(o: s Aqcd/toq) and 0(a s aAQCD), replacing 
Oy* in eq. ([!]) by an alternative set of operators 
O'y — hTyq hTyq is equivalent to eq. (Q) itself 
up to the higher order corrections. 

In the following, we explain our method by tak- 
ing Bp as an example. First we define the "lattice 
_B-parameter" by 
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where A^^joj^q. B' Y a is also defined in 



similar manner but with O'y and A' 1 ^ =h r y^^q. 
The Bb is, then, calculated in the following four 
forms, 

1: B B = {B^ + £ y z LX B l f) /(l + z A f 

2: B B = + J2y z L.Y-A 2 By t 

3: B B = (b^ + £y z L}Y B'f) /(l + z A ) 2 
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where za is the one-loop matching coefficient of 
axial vector current and Zjj y—a 3 = z l,y except 
for zj / j,__42 = zl_l — 2za- The difference between 



the method 1 and 2 (or between the method 3 and 
4) is of 0{a 2 s ) : since in the former expression the 
factor 1/(1 + za) 2 , which arises as the renormal- 
ization factor of Aq in the denominator, is kept 
untouched and in the latter it is expanded to leave 
only the linear term in a s . On the other hand, 
the difference between the method 1 and 3 (or 
between the method 2 and 4) is of 0(a s p/m,Q) or 
0(al/am,Q), since the former uses Oy* and the 
latter O'y ■ Note that all these definitions of Bb 
have equal accuracy up to higher order contribu- 
tions as just discussed. A typical size of these 
uncertainties may be estimated from scatter of 
the results. For the coupling constant, which is 
included in Zx,y, <xv(q*) with q* = 2/ a is used 
throughout this work. 

In the case of B$, the constant 8/3 and axial 
vector current Aff* in the denominator of eq. @), 
za and zl y^a 2 are replaced by 5/3 and pseudo- 
scalar current P lat = b'fsq, zp and zs.y-p 2 , re- 
spectively. 

3. Results 

In the following, we show the results of 
BB d (rrib), Bs s /BB d and i?5 s (mb) obtained with 
the four methods. Figure [l] shows the 1/Mp de- 
pendence of BB d {nib) obtained at /3=5.9, where 
Mp is the pseudo-scalar heavy-light meson mass. 
The renormalization scale of the continuum op- 
erator is set to ^i=TOb=4.8 GeV. First we find 
that four plots in Fig. [I] show various 1/Mp de- 
pendences and there is no clear tendency. This 
means that Bp depends on 1/Mp only weakly, at 
least, in the region of rriQ > 2to{,/3. Next one 
can see that the difference between the method 1 
and 2 (or 3 and 4) becomes larger as one goes to 
the static limit and dominate the scatter of the 
results near the static limit. This might be inter- 
preted as a typical size of the 0(a 2 s ) uncertainty 
as discussed before. We also find a deviation at 
the lightest Mp, which is dominated by the dif- 
ference between the method 1 and 3 (or method 
2 and 4). At this lightest Mp, the bare heavy 
quark mass in lattice unit is not much larger than 
one (amQ=1.3), so we might infer that this dif- 
ference indicates a typical size of the 0(a 2 /arriQ) 
error rather than the 0{a s p/mQ). Interpolating 
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Figure 1. 1/Mp dependence of B Bd (mi,) at = 
5.9. 

all the data to the physical B meson mass, typical 
sizes of the uncertainties are estimated as 



0(a 2 s ) ~ 5% 
0(a 2 s /amQ) + 0(a s p/rriQ) 



4% 



(4) 



Repeating the same analysis at three lattice spac- 
ings, we obtain the lattice spacing dependence of 
BB d , which is plotted in Fig. [|. From this fig- 
ure we can see that the three data are consistent 
with each other within the magnitude of scatter 
of the results. This suggests that the scaling vio- 
lation for B Bd is not too large. We take our best 
estimate from the central value of the results at 
/3=5.9. The remaining uncertainty is estimated 
by the naive order counting, which gives 



0(a 2 s ) ~ 7%, 0(a sP /m b ) ~ 
0(a 2 Jam b ) ~ 4%, 0(a 2 p 2 ) 



2%, 

~ 3%, 0{a s ap) 



5%. 



Adding them in quadrature we obtain 10% as our 
estimate for the systematic error, which is more 
conservative than eq. (Q), that is, just taking scat- 
ter of the results as systematic uncertainty. Our 
preliminary result in the quenched approximation 
is i3_B d (wb)=0.84(2)(8), where the errors are the 
statistical and systematic ones, respectively. 

Repeating the similar analysis, we obtain the 
lattice spacing dependences of B B JB Bd shown 
in Fig. |^. For this quantity, we observe neither 
method dependence nor large scaling violation. 
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Figure 2. BB d {mb) as a function of lattice spac- 
ing. 

We take our best estimate from /3=5.9 lattice as 
before, and obtain B B jB Bd = 1.017(10) (t£) by 
only considering the uncertainty in n s . 

Figure ^ shows the lattice spacing dependence 
of Bs s ■ The data plotted in this figure does not 
show a large scaling violation within the system- 
atic uncertainty indicated by the scatter of the 
results. Following the same analysis as in the 
case of B Bdl we obtain B s , =0.87(1) (9) (to) ™ 
the quenched approximation, where the errors are 
statistical, systematic and coming from the uncer- 
tainty of k s , respectively. 

4. Applications 

Let us now consider a few applications of our 
results. The first one is the mass difference in the 
neutral B q meson systems, which is given by the 
CKM matrix elements (|V^Vt g |) and f B B Bq up 
to a known factor as 
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(5) 



Assuming the additional 10% uncertainty due to 
the quenched approximation to above result for 
B Bd and the recent world average of the de- 
cay constant from n/=2 simulations /s d =210(30) 
MeV, we obtain |V^V W |=0.0074(12). This result 
can be used to estimate the frequency of the B° s - 
B° s mixing. If we take |V ts HV cfc |=0.0404(18) and 
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Figure 3. Bb s ) Bb a as a function of lattice spac- 
ing. 

/s s //s (i =1.16(4) and assume the uncertainty due 
to the quenched approximation for BB s /BB d to 
be 10% again, we obtain AM Ba = 20.0(4.7) ps" 1 . 
The current experimental bound is AMb 3 > 
14.9 ps- 1 [@. 

The second application is the width differ- 
ence of B s mesons. The width difference is 
theoretically calculated using the 1/rrib expan- 
sion. Beneke et <xl. have completed the for- 
mula with the next-to-leading order correction 
in Ref. ||. Following their analysis but with 
/ B ,=245(30)MeV H, the formula is given by 



AT 

r 



h 



-X 



0.008B Bs (m b ) 



245(30)MeV J 
0.227(17)S Ss (m 6 ) - 0.086(17) . (6) 



Using the results for i3-parameters but with the 
additional 10% error due to the quenched approx- 
imation, we obtain (Ar/I% = 0.119(29) (32) (17), 
where errors come from the uncertainties of fs s , 
Bs s and the last term in eq. @, which shows 
0(l/mb) contribution, respectively. The current 
experimental bound is (AT/T) S = 0.16(±J|) |]. 
It would be quite interesting to compare our re- 
sult with the coming Tevatron Run II data, which 
are expected to improve the current experimental 
results significantly soon. 
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Figure 4. Bg s (jnj,) as a function of lattice spacing. 
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